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Abstract 

Two-dimensional Heisenberg model with anisotropic couplings in the x and y 
directions ( J x ^ J y ) is considered. The model is first solved in the Schwinger- 
boson mean-field approximation. Then the solution is Gutzwiller projected 
to satisfy the local constraint that there is only one boson at each site. The 
energy and spin-spin correlation of the obtained wavefunction are calculated 
for systems with up to 20 x 20 sites by means of the variational Monte Carlo 
simulation. It is shown that the antiferromagnetic long-range order remains 

down to the one-dimensional limit. 
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I. INTRODUCTION 



Stimulated by the discovery of cuprate superconductors, the ground state of the two- 
dimensional S =1/2 quantum antiferromagnetic Heisenberg model has been extensively in- 
vestigated by various kinds of approaches. Although the model has never been solved exactly, 
the properties of its ground state have been carefully investigated by numerical and analyti- 
cal calculations. It has been found that ground-state energy per site is (—0.6696 ± 0.0004) J 
by series expansions!, and (—0.6692 ± 0.0002) J by Green's function Monte Carlo method!. 
As for the existence of the long-range order, it is believed that it has an antiferromagnetic 
long-range order. The staggered magnetization per site is estimated to be 0.308 ± 0.008 by 
series expansions^! and 0.31 ± 0.01 by Green's function Monte Carlo method.! 

On the other hand, one-dimensional S = 1/2 quantum antiferromagnetic Heisenberg 
model is solved exactly by means of the Bethe ansatz. The exact ground-state energy per 
site is (In 2 — 0.25) J = — 0.4431J and the excitation is gapless.! Since the system is one- 
dimensional, it does not have a long-range order. 

From these results, we expect an order- disorder phase transition in an anisotropic Heisen- 
berg model which interpolates the one- and two-dimensional systems; 

W = J X Sj • Sj + £ + Jy Sj • Sj_|_y, (1) 

i i 

where i + x (i + y) means the site next to the z-th site in the x(y) direction, and J x = J 
and J y = aJ are the coupling constants in the x- and y-direction, respectively. Since 
the ground state of this Hamiltonian has no long-range order at a = (one-dimensional 
case) and it has the antiferromagnetic long-range order at a = 1 (two-dimensional isotropic 
case), we can expect a quantum phase transition as a is changed between and 1. It had 
been believed that this phase boundary is at a = 0. When a is infinitesimally small, the 
system can be regarded as a collection of weakly-coupled one-dimensional chains. A simple 
mean-field approximation indicates that infinitesimally small a makes the system ordered.!! 
However, very recently Parola, Sorella and Zhongiil asserted that they found an evidence 
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for an order- disorder transition at finite value of a. They investigated the magnetization 
by spin wave theory and Lanczos method. The spin wave theory predicted a breakdown 
of the antiferromagnetic long-range order at about a c ~ 0.03(0.07) at first order (second 
order) in 1/25 expansions. Below this critical anisotropy, the staggered state disappeared 
and disordered state was realized. They also found that the spin wave theory generalized 
to finite systems agreed with the results of the Lanczos method for a > 0.1. However, this 
agreement did not persist for a < 0.1. The result of the Lanczos method showed finite 
magnetization at a = 0. Although they concluded the order- disorder phase transition at 
finite a from their results, it is far from conclusive. 

In this paper we investigate this order-disorder transition using a different method. Here, 
we consider a resonating valence bond (RVB) stateffi and adopt the Schwinger-boson de- 
scription of the spin operators. For the isotropic two-dimensional case (a = 1), wave 
functions based on the RVB singlet have been studied by many authors. Liang, Dougot 
and AndersonS made a variational RVB wave function with long-range bonds and found 
that the optimal energy per site is —0.6688 J and the ground state has an antiferromagnetic 
long-range order. The staggered magnetization is 0.225 per site. The variational energy is 
very close to the best estimated value.0 Furthermore, recently Chen and XiuQ made a trial 
wave function by the Gutzwiller projection of the state obtained in the Schwinger-boson 
mean-field theory. Though the wave function has no variational parameter, they found that 
it can describe the isotropic ground state quite accurately. The energy per site is — 0.6688J 
which is 99.9 % of the best estimation and the staggered magnetization is also very close 
to the best estimated value. In this paper, we extend the wave function of Chen and Xiu 
into the anisotropic Hamiltonian for < a < 1 to investigate the order-disorder transition. 
To this end, we first solve the model in the mean-field approximation. Next, using the 
mean-field solution, we make variational wave functions without doubly-occupied sites and 
calculate the energy and staggered magnetization to investigate the phase transition. 

In section [D], we solve the Hamiltonian in the Schwinger-boson mean-field theory. Using 
this mean-field solution, we make RVB wave functions without doubly occupied sites in sec- 



tion |ITT| . In section |V[ numerical results for energies and spin-spin correlations as a function 
of anisotropy are studied by means of Monte Carlo simulations. In the last section [V], we 
discuss accuracy of our wave function and the boundary of the order- disorder transition. 



II. MEAN FIELD SOLUTION 



In this section we solve the Hamiltonian in the Schwinger-boson mean-field 

T and 



approximation.l^Jfi^M First we introduce two kinds of bose operators Sj t and Sj i to express 



the spin operators, 



S t = s ^i,1 s i,b and S i = \( s h s ^ ~ s l,i s u)- ( 2 ) 



The commutation relations of the spin operators Si are satisfied in this replacement. We 
impose a constraint, 

s h s i,r + 4,i s i,i = !> (3) 

in order to guarantee S = 1/2. Then the Hamiltonian is rewritten as follows: 

i z ^ i 

+A t H H S^i^Si^ . (4) 

Here z = x or y, i+z represents a site next to the site i in the ^-direction, n,i = sj+s^f+sj rSv[_, 
and /i is a chemical potential to be adjusted to satisfy Eq.([|) on the average. 

To solve this Hamiltonian in a mean-field approximation we introduce the following 
mean-field parameters which give the amplitude of the nearest-neighbor singlet pairs and an 
averaged occupation number, 

&z = ^{s u Si +S)1 - SijSi+^i), (5) 

n* = (4,a S i,a) = ( 6 ) 



After decoupling the Hamiltonian, we rewrite the operator by its Fourier transform: 
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where N is the total number of lattice sites, and k are two-dimensional wave vectors, 
mean-field Hamiltonian is then written as 

Ti MF = [ A ( s k,T s k, T + s f -ka s -k,i) + 7k4,T st -k,| + 7 k S-k,|Sk,T] 

k 

+ 2A^[J|A X | 2 + Ja\A y \ 2 ] + ±N{1 + a) J, 

where 

\ = H- ^(1 +a)J, 
7k = 2i( JA X sin k x + JaA y sin k y ). 

This Hamiltonian is diagonalized using the Bogoliubov transformation: 

Sk,t = u k a k - v k (3l, 

s-k,i = ~v k al + w k /3 k , 

where 

«k = -^[^ + l]^xp(^ k ), 



and 



^ = -^[^-l]^exp(^ k ). 



#k is the phase of 7^, 7^ = |7fc| exp(i0 k ), and 



k 



E k = v/A 2 - |7 
After this transformation the Hamiltonian becomes 
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T-Cmf = E k( a i a k + PkPk) + const., (16) 

k 

The ground state \G) is defined as the vacuum of the bose operators and (5y : a k \G) = 
(3)i\G) = 0. For a finite size system, A, A^, and A y are determined from the following self- 
consistent equations, 

(G|4 T ^ )T + 4^4^) = 1, (17) 



and 



which read 



(G|si4Si+f )T - Sij8 i+ t t i\G) = 2A Z , (18) 



2 4?i < 19 > 



and 



. 1 x — *s , Jr^^v sin /u-r. — I - J<i,^^fii sin /u7< 

A * = 77 £ sin fc*"^ F ( 2 °) 

The solution depends on the size of the system N, so does Uk and f^. We will use the 
solutions for finite size systems in the following sections. 

When N is finite E^ never becomes zero. However, in the limit of N — > oo it is possible 
that Ek vanishes at k = K.± = ±(tt/2, tt/2). In such a case it is known that we need to 
introduce the Bose condensate nJHHll, and Eqs. (|19|) and ( p0|) are rewritten as 

2 =(2W-W->4 +2nfl ' (2l) 

a, = JL_ £ dkl £ ik< sm ^A-u^*t, + „ b (22) 

Here becomes finite only when the spectrum is gapless, i.e. when A = 2(J X A X + J y A y ). 
The finite value of ub means the existence of the antiferromagnetic long-range order. The 



staggered magnetization is y 3/2 

Before discussing the mean-field solution for an arbitrary anisotropy, it is worth while 
to see what happens in the one-dimensional limit where a = 0. It is easy to see that the 



spectrum cannot be gapless, and thus tib = 0. In this case the self-consistent equations 
are expressed by elliptic integrals. By introducing e = 2JA X /X, we can express Eq.(|2T|) and 
Eq.([22p as 

2 = J- r &a = 2 K(e); (23) 

and 

A 1 r „ sin 2 ^ 



= 4"A E ( e )- (24) 

Here K(e) and E(e) are the complete elliptic integrals. By solving these equations we obtain 
A = 1.379962J, and A x = 0.6792397. The excitation energy E k has a gap: E k > 0.242547J. 

For intermediate value of a we must solve the self-consistent equations numerically. The 
results are shown in Fig.l, where the RVB order parameters A x and A y , the bose condensate 
tie, and the energy gap E g are plotted. When a is decreased from unity, both and 
decreases. The bose condensate vanishes first at a = 0.1356. This is the point where the 
antiferromagnetic long-range order vanishes. Below this value of a, the energy gap in the spin 
excitation develops quite rapidly and A y is suppressed. A y finally vanishes at a = 0.1286. 
Below a = 0.1286 the system reduces to a collection of independent chains running parallel 
to the x-direction. In this region the order parameter A x and the gap E g become constant. 

III. RVB WAVE FUNCTIONS 

The ground state wave function obtained in the mean-field theory is expressed as, 

|G) = nex P f-^4,A,l)l°)' (25) 

where |0) is the vacuum of the Schwinger-bosons. By the Fourier transformation for sj.^ , 
I' 



, , we get a real space representation for this ground state; 
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|G) = expEa M 4 T 4|]|0), (26) 

„ = --u- 



a 



H^expJ-ik.ry), (27) 



where r^- = r^ — r^. It is evident that the local constraint, Eq.@, is not satisfied in this wave 
function. This is why the energy of the mean-field ground-state is too low .111 We remove 
this difficulty by projecting the wave function to a space where every site is singly occupied. 
Namely, we perform the Gutzwiller projection, 



|G> = ^(E^Wi)'l°>- ( 28 ) 

The operater P in the right hand side is the Gutzwiller projection operater. From now on, 
we consider this wave function. Since a^j = — a^j, the ground state of Eq.(^Sj) is nothing 
but an RVB state which includes long-range bonds with weights a^j. 

Although it would be possible to regard every djj as a variational parameter, we here 
restrict a^j to be those given in Eq . (p7|) . In the isotropic case (J x = J y ), this restriction is 
justified by the result itself: Chen and Xiu0 have shown that this choice of Ojj gives excellent 
results for the ground-state energy and the staggered magnetization. To see if Eq.(j27]) gives 
similarly good results for anisotropic case or not is one of our aims in the present paper. It 
should be noticed that for the anisotropic system we still have one variational parameter. 
Namely, and Vk depend on the anisotropy parameter a = J y / J x . It is not evident that for 
a given value of a in the Hamiltonian, the same value of a in the variational wave function 
gives the best result. Thus we consider a in the mean-field equations as a variational 
parameter which we denote a p . 

IV. NUMERICAL RESULTS 

In this section, we show numerical results of the variational energy, spin-spin correlation, 
and staggered magnetization for the anisotropic Heisenberg Hamiltonian. We perform Monte 
Carlo simulations for lattices with various number of sites up to 20 x 20 for the energy 
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and the spin-spin correlation. All the numerical calculations are performed with periodic 
boundary conditions. For each system size we solve the self-consistent equations (|1^) and 
(p0|), and calculate djj to be used to construct the wave function at that system size. Instead 
we could have used aij for an infinite size system. However, we did not take this approach 
since the energy is higher and size-scaling does not coincide with the prediction of the spin 
wave theory in this case. 

A. Ground state 

The energy per site of the anisotropic Heisenberg model is 

E = J x (e x ) + J v (e y ), (29) 

where 

(^H^Em-s^G). (30) 

i 

Here (e z ), z = x or y, is an expectation value of the nearest-neighbor spin-spin correlation in 
the x- or y-direction. These (e x ) and (e^) depend on the system size N and the variational 
parameter a p , but independent of a. Thus we first estimate the thermodynamic limit of (e^) 
and (e y ) for several values of a p . The size dependence is studied and we find the following 
size scaling 

(e z (L)) = (e z (oo)) + X Z L- 3 + ■ ■ ■ , (31) 

where A^, is a constant and L is the linear dimension, L 2 = N. The size-scaling coincides 
with the spin wave theory for a square lattice. Then we obtain a p dependence of the energy 
E at a fixed a for the infinite size system. Fig. 2 shows such dependence at a = 0.6 for 
example. The error bars result both from Monte Carlo statistical errors and from a fitting 
error of the size-scaling. It can be seen that the optimal energy is given when a p = 0.680. In 
a similar way, we obtain the optimal energy for various values of a, which is shown in Fig. 3. 
When a is equal to zero, the system becomes the one-dimensional S = | antiferromagnetic 
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Heisenberg chain. In this case, the exact energy per site is given in the Bethe ansatz 
solution which is (In 2 — 0.25) J = — 0.4431 J.! Our variational energy per site in this case 
is (— 0.4337 ± 0.0030) J which is 97.9% of the exact result. This shows that our RVB wave 
function can also describe fairly well the anisotropic Heisenberg Hamiltonian as well as the 
isotropic one. 



B. Staggered magnetization 

There are several methods to estimate the staggered magnetization. In most numerical 
calculations, it is estimated from the spin-spin correlation. Here, we use the staggered 
spin-spin correlation at the longest possible distance for each direction,0 

Ml{L)=C z {^\ (32) 
(-l) n 

^) = [^H£<|S«-S^|>, (33) 

i 

where z = x or y and % + nz is the n-th site from the site i in the x- or y-direction. We 
calculate M Z (L) for each lattice size and extrapolate to the thermodynamic limit using the 
finite size scaling: 

M Z (L) = M z (oo) + fiL" 1 + ■ ■ ■ . (34) 

where \i is a constant. This scaling agrees with the prediction of the spin wave theory and 
arguments given by Huse et al.@. The obtained M 2 (oo) is the staggered magnetization. 

Since the variational function is a Gutzwiller projected mean-field solution, the properties 
of the latter wave function is inherited by it. The mean-field wave function which shows the 
long-range order gives a finite staggered magnetization even after the Gutzwiller projection, 
while that which shows no long-range order gives vanishing staggered magnetization after 
the projection. We find that, for any anisotropy, the optimized parameter a p is larger 
than 0.1356 for any a and, as a result, our variational state always has a long-range order. 
Furthermore we find that the magnetizations in the x- and ^-directions coincide, which is 
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consistent with the mean-field solution. The a dependence of the staggered magnetization 
is shown in Fig. 4. The staggered magnetization and hence the antiferromagnetic long-range 
order remains down to a = 0. 

V. DISCUSSIONS 

In this paper we first solved the Hamiltonian by the Schwinger-boson mean-field theory, 
and then the solution is Gutzwiller projected to obtain the variational wave function, which 
we used to investigate the order-disorder transition. Our Gutzwiller-projected wave function 
is very suitable for studying the order- disorder transition. Since it is rewritten as a RVB 
state, it can represent a disordered spin state as well as an ordered state. As shown by Liang, 
Dougot and AndersonEl, when the weight of the singlet RVB bond, a^-, has a power-law 
behavior, the RVB state has a long-range order. On the other hand, when cjjj is short- 
ranged, the RVB state becomes a disordered state with a spin gap and it has a short range 
spin-spin correlation. 

By solving the Schwinger-boson mean-field theory for anisotropic Heisenberg model, we 
construct anisotropic djj which is suitable for the anisotropic case. Furthermore, we can 
compare the variational energy for the ordered state (a p > 0.1356), and the disordered state 
(qj p < 0.1356) by regarding Q^p clS cl variational parameter. We have found that the optimized 
parameter a p is always larger than a and moreover larger than a c = 0.1356. (Even at a — > 0, 
we get a p = 0.205.) This means that the ground state is stabilized if it develops a long- 
range order for all parameters, < a < 1. Although both the spin wave theory discussed 
by Parola et al.i and the Schwinger-boson mean-field theory predict a disordered phase for 
small a, our result supports the existence of a long-range order even for small a. 

We believe the phase transition obtained in the spin wave theory and in section |I| at 
finite value of a is an artifact of the poorness of the mean-field solutions. It is known 
that in the one-dimensional limit even though there is no long-range order, the excitation 
spectrum is gapless and the spin-spin correlation decays not exponentially but algebraically. 



11 



This algebraic decay lies in the base of our belief that an infinitesimal coupling between one- 
dimensional chains causes a long-range order in the entire system although the magnitude of 
the staggered magnetization being infinitesimally small. On the other hand our mean-field 
solution at J y = gives a gapfull excitation spectrum, and the spin-spin correlation shows an 
exponential decay. This qualitative difference between the mean-field and the exact solution 
does not persist in our variational wave function, since we use a p > 0.205. We remark here 
that the gapfull excitation and the exponential decay of the spin-spin correlation are the 
properties of the integer spin antiferromagnetic Heisenberg chain, and the order-disorder 
transition in the presence of the interchain coupling has been successfully discussed by 
Azzouz and DougotEl by the Schwinger-boson mean-field theory. 

Our wave function is much improved from the mean-field solution by the Gutzwiller 
projection and by the optimum choice of the variational parameter a p . However, the energy 
and the magnetization in the limit of a — > clearly indicates that there is ample room 
for improvements of our wave function. Our results show discontinuity in the staggered 
magnetization at a = 0: with infinitesimal a the magnetization jumps from zero to 0.1227± 
0.023. This will not be correct. When the ordered state persists down to the limit of a = 0, 
the most natural behavior would be a continuous growth of the magnetization from zero as 
a increases; power dependence on a being quite plausible. We expect improvements of the 
ground-state energy at a < 1 will decrease the value of magnetization. At present we do 
not know if such an improvement will cause the magnetization vanish even at finite value of 
a or not. However, our results show strong tendency of the existence of a long-range order 
even for infinitesimally small a. 

In conclusion we have investigated the two-dimensional anisotropic Heisenberg model by 
the Schwinger-boson Gutzwiller projection method. Our results show existence of the long- 
range order in the whole range of the anisotropic parameter a. This is in disagreement with 
the work by Parola et al.i We think the order-disorder transition at finite value of a is still 
an open question. More elaborate work will be necessary to draw a definite conclusion for it. 
We can and will try to improve the wave function by treating a^j as variational parameters. 
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The results will be reported in the near future. 
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FIGURE CAPTIONS 



Fig. 1. Mean- field values of order parameters A x ,A y , bose-condensate n B and energy 
gap E g as a function of a. 

Fig. 2. Variational energy per site as a function of a p for a = 0.6. The optimal energy 
is given at a p = 0.680. 

Fig. 3. Energy per site as a function of a. When a — 1, isotropic case, the energy is 
(-0.6688 ±0.0008) J, and when a = 0, one-dimensional case, it is (-0.4337 ± 0.0005) J. The 
exact value at J y = is shown by an arrow. 

Fig. 4. Staggered magnetization M = M x (oo) = M y (oo) as a function of a. 
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